clear all
set more off
set scheme s2color
cd ""

use "sample.dta", clear

bysort prov_code: keep if _n == 1

#delimit ;
histogram prov_age, width(0.5) density
          xtitle("Age (years)") ytitle("Density")
          xlabel(, labsize(small)) color(orange%60)
          ylabel(, grid tlength(0.8) glwidth(vthin) glpattern(dash) nogextend labsize(small))
          graphregion(fcolor(white) color(white) icolor(white) ilstyle(p1box)
          margin(small)) plotregion(margin(vsmall));
#delimit cr

replace prov_tenure = prov_tenure - 5 if prov_tenure > 5
#delimit ;
histogram prov_tenure, width(0.5) density
          xtitle("Tenure length (years)") ytitle("Density")
          xlabel(, labsize(small)) color(orange%60)
          ylabel(, grid tlength(0.8) glwidth(vthin) glpattern(dash) nogextend labsize(small))
          graphregion(fcolor(white) color(white) icolor(white) ilstyle(p1box)
          margin(small)) plotregion(margin(vsmall));
#delimit cr
